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The possible formation of oxides or thin oxide films (surface oxides) on late transition metal 
surfaces is recently being recognized as an essential ingredient when aiming to understand catalytic 
oxidation reactions under technologically relevant gas phase conditions. Using the CO oxidation at 
Pd(100) as example, we investigate the composition and structure of this model catalyst surface over 
a wide range of (T,p)-conditions within a multiscale modeling approach where density- functional 
theory is linked to thermodynamics. The results show that under the catalytically most relevant 
gas phase conditions a thin surface oxide is the most stable "phase" and that the system is actually 
very close to a transition between this oxidic state and a reduced state in form of a CO covered 
Pd(100) surface. 

PACS numbers: 82.65.+r, 68.43.Bc, 68.43.De, 68.35.Md 



I. INTRODUCTION 

Catalytic oxidation using transition metals (TM) as 
the active material is an important technological process, 
of which we still have only limited microscopic insight. 
Much valuable microscopic information that has been ob- 
tained under ultra-high vacuum (UHV) conditions is not 
directly transferable to the conditions of practical cataly- 
sis (with pressures of several bars and elevated tempera- 
tures). And in-situ techniques that would operate under 
these conditions are still struggling to achieve atomic- 
scale surface sensitive information. A key factor that hin- 
ders the direct transfer of the UHV insight is that under 
reaction conditions the entire structure and composition 
of the catalyst surface might be changed (often called the 
"materials gap" ) . And if a new material is created during 
the induction period of the catalytic process, it will obvi- 
ously exhibit a very different chemical activity; it will be 
quantitatively and even qualitatively different. Under the 
oxygen-rich gas phase conditions and elevated tempera- 
tures of oxidation catalysis the surface of the transition 
metal catalyst might e.g. be oxidized, so that the actual 
catalytically active state is not the pristine metal, but 
rather the formed oxide. One example for this is the CO 
oxidation reaction over ruthenium. Here it was found, 
that under UHV conditions the Ru(0001) model catalyst 
shows almost no catalytic activity, whereas under high 
oxygen pressures its catalytic activity exceeds even the 
one of the frequently used palladium and platinum cata- 
lysts 1 . This increase in the catalytic activity was traced 
back to the formation of a Ru02 film at the surface^. 
For the ruthenium catalyst a clear distinction could thus 
be made between the different states, namely the weakly 
catalytically active metal under low oxygen pressures and 
the highly active oxide film under high oxygen pressures. 

For the here discussed model system, the CO oxidation 
at Pd(100), recent in-situ experimental measurements 
also indicate that an oxidic structure at the Pd(100) sur- 



face might be formed under the gas phase conditions ap- 
plied in industrial oxidation catalysis. Using surface x- 
ray diffraction (SXRD) measurements a (T, p)-diagram 
showing the detected structures in a pure oxygen atmo- 
sphere over a temperature range of T = 600 — 1000 K 
and pressures of po 2 = 10 -9 — 1 atm could be inferred 4 . 
A thin surface oxide structure was measured over an ex- 
tended (T,p)-range, suggesting that this structure might 
also appear under catalytic reaction conditions. In re- 
actor scanning tunneling microscopy (STM) experiments 
by Hendriksen et al 5 ' 6 the structure of the Pd(100) sur- 
face could directly be monitored during the catalytic ox- 
idation of CO. Here, the surface was exposed to both 
oxygen and CO at a total pressure of p tot « 1 atm and a 
temperature of T w 400 K. The partial pressures of the 
reactant gases CO and O2, as well as the reaction prod- 
uct CO2, were measured simultaneously with the STM 
images. Depending on the partial pressure of O2 and CO 
the authors observed a change in the reaction rate, which 
was accompanied by a significant change in the morphol- 
ogy of the surface. This morphology change (albeit not 
atomically resolved) was interpreted as a change from the 
adsorbate covered Pd(100) surface to an oxidic state. 

To address this problem from a theoretical point of 
view we employ a multiscale modeling approach where we 
use density-functional theory (DFT) to describe the sys- 
tem on an electronic (microscopic) level. By linking these 
results to thermodynamics^*^ it becomes possible to 
address much larger system sizes and to compare the sta- 
bility of different surface structures in contact with the 
surrounding gas phas o n i 12 i 13 i 14 i 15 i 16 . We use in partic- 
ular a constrained equilibrium approac h 15 ! 16 , where the 
surface is considered to be in full thermodynamic equi- 
librium with two separate gas phase reservoirs of O2 and 
CO, and not allowing that O2 and CO can react. The 
key result obtained is a surface "phase diagram" , which 
provides first insight into possible surface structures over 
a wide range of temperature and pressure conditions of 
the O2 and CO gas phase. Focusing on gas phase con- 
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ditions relevant for oxidation catalysis we find the sys- 
tem on the verge of either stabilizing a thin surface oxide 
structure or a CO covered Pd(100) surface. It is thus 
well possible and even likely that the surface oxide struc- 
ture contributes to the active state of the Pd(100) model 
catalyst under reaction conditions. However, to verify 
this, the kinetic effects of the on-going catalytic CO2 for- 
mation need to be considered. This is done in a second 
step of our hierarchical approach by performing kinetic 
Monte Carlo (kMC) simulations on the here identified 
surface structures, on which we report in a consecutive 
paper. 



II. THEORY 

A. Gibbs free energy of adsorption 

To determine the stability of a surface in contact with 
a gas phase reservoir we use the surface free energy 7. For 
a multi-component system in equilibrium with atomic 
reservoirs (e.g. a surrounding gas or liquid phase environ- 
ment, or a macroscopic bulk phase) a general expression 
for the surface free energy is given by 



l(T,{Pi}) = \ 



(1) 



Here, G surf is the Gibbs free energy of the solid includ- 
ing the surface, A is the surface area and fii(T,pi) are the 
chemical potentials of the various species i in the system. 
To compare the stability of different adsorption struc- 
tures of oxygen and CO on the Pd(100) surface depending 
on the surrounding gas phase conditions one has to eval- 
uate the Gibbs free energy of adsorption, AG ads For 
this, the stability of the different adsorption structures is 
compared with respect to the clean metal surface. For a 
Pd(100) surface in a constrained equilibrium with an O2 
and CO gas phase AG ads is thus given by 
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where Gq co@Pd(ioo) * s ^ ne Gibbs free energy of the 
Pd(100) surface with No adsorbed O atoms and Nqo 
adsorbed CO molecules, and Gp U d r ^ f 100 ^ is the Gibbs free 
energy of the clean palladium surface. If the number 
of palladium atoms per surface area A in the adsorp- 
tion structure, iVpd, and in the clean surface, A/p d , are 
not equal (e.g. due to a surface reconstruction), i.e. 
ATVpd = TVpd — Np d 7^ 0, then the excess/deficiency 



atoms are taken from/put into a bulk reservoir, repre- 
sented by the chemical potential of the bulk solid phase, 
/ip d lk . The chemical potentials of the two gas phase reser- 
voirs of oxygen, /io 27 an d CO, /ico, have been separated 
into a total energy contribution and the remaining part 
containing all the temperature and pressure dependent 
terms, i.e. /i = l/2/xo 2 = 1/2^* + A/i (T,p), set- 
ting A/i (T,p) = l/2A/i 02 (T,p), and /i CO = £g£ + 
A/ico (T,p)iil£. In the last line of Eq. © the differ- 
ence of the Gibbs free energies of the clean and adsorbate 
covered surfaces, and of the bulk system, has been ap- 
proximated by the DFT total energies E^^. A detailed 
discussion of the thereby neglected contributions to the 
Gibbs free energy of adsorption and their influence on 
the presented results will be given below. 

In the constrained equilibrium approach leading to 
Eq. ([2]) the two gas phases are assumed to be non- 
interactingi^. The Gibbs free energy of adsorption AG ads 
depends then linearly on the chemical potentials of both 
the oxygen and CO gas phase. The slope is respectively 
determined by the surface coverage of oxygen and/or CO. 
This is most apparent in the last line of Eq. ([2]), where 
we have introduced the average binding energy of O and 
CO as 
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Since the most stable structure will be the one with 
the lowest surface free energy, an adsorbate structure 
will be stable with respect to the clean surface, if 
7o,co@Pd(ioo) < 7Pd(ioo) ? i- e - if AG ads > 0. If two struc- 
tures contain an equivalent amount of oxygen and CO 
they will also show the same dependence on A/io and 
A/ico- The ratio of how these two structures contribute 
to the stable structures is governed by the law of mass ac- 
tion, which enables us to directly exclude a less favorable 
structure, if its binding energy, ^oco@Pd(ioo) > differs by 
much more than k^T from the one of the most favorable 
structure. 

For the adsorption of O/CO on the Pd(100) surface the 
average binding energy defined in Eq. (j3]) is equivalent to 
the commonly used binding energy, which is often given 
per O atom respectively CO molecule, e.g. in the case of 
oxygen as 
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In the case of additional adsorption of O/CO on the re- 
constructed y/h surface oxide (which will be introduced 
below) the binding energy of the atom/molecule with re- 
spect to this y/b structure is then given by (here exem- 
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where E^ is the total energy of the y/E surface oxide. 

We notejjhat E hind is in this case different to the above 
defined E hind in Eq. ©, since the latter also contains 
the energetic changes due to the formation of the y/h 
surface oxide, i.e. it is referenced with respect to the 
clean Pd(100) surface. 



B. Bulk oxide stability 



TABLE I: Computed heat of formation Aif| dO (0,0) and 
binding energies of O2, CO and CO2 (see text). The quoted 
experimental values are extrapolated to T — K and the zero 
point vibration energy (ZPE) has been removed. All values 
are in eV. 





PBE 


RPBE 


LDA 


Exp.^ 


ZPE 


Afl* dO (0,0) 


-0.87 


-0.62 


-1.42 


-0.97 




7-ibind 

7-ibind 
^CO 
TTibind 

^co 2 


-6.22 
-11.65 
-17.99 


-5.75 
-11.20 
-17.09 


-7.56 
-12.93 
-20.43 


-5.30 
-11.33 
-16.98 


(0.10) 
(0.13) 
(0.27) 


AE mo1 


-3.24 


-3.02 


-3.72 


-3.00 





When considering a metal surface in contact with an 
O2 and CO gas phase complete conversion into an ex- 
tended bulk oxide is also a thermodynamic possibility. 
The stability of the corresponding bulk oxide has there- 
fore to be evaluated with respect to the two gas phase 
components and compared to the various stable struc- 
tures forming at the metal surface 18 . A sufficient oxygen 
content in the gas phase will lead to the formation of bulk 
palladium oxide, PdO, whereas a sufficient CO content 
will favor the decomposition of PdO into CO2 and Pd 
metal. In a pure oxygen gas phase the thermodynamic 
stability of the bulk oxide is given by 

..bulk ^ .bulk _j_ (n\ 
MPdO < MPd + Mo , (o) 

whereas in a pure CO environment the stability condition 
for the oxide is 

Mpdo + Mco < Mpd lk + Mcos • (7) 

If we approximate the chemical potential of the free CO2 
molecule by its total energy only 18 , Eq. ([6]) and Eq. 
can be combined yielding the stability criterion for PdO 
in an oxygen and CO containing gas phase 

A/ico-A/i O <-2A^ dO (T = 0K) + A^ mo1 . (8) 

Here, AH* d0 « E$% Q - - 1/2E$* is the heat of for- 
mation of PdO at T = K, and AE mo1 = Egg* - E$£ d - 
l/2£ , Q 1 ^ ld is the difference in binding energies between the 
three gas phase molecules. Table U compiles these bind- 
ing energies, as well as the heat of formation of PdO, as 
obtained by previously published highly converged DFT 
calculations^ and using the computational methodology 
described below. In addition the experimental values are 
quoted. These values have been extrapolated to T = K 
and the zero point vibration energy (ZPE, also listed in 
Table [J) has been removed. 

Substituting the chemical potential of CO2, /ico 2 ? by 
only the total energy term, Eqq , in Eq. is a rather 
crude approximation. The most important contributions 
that have been neglected arise from the vibrational and 



translational free energy. The assumption that the CO2 
formed at the surface is readily transported away, moti- 
vates us to disregard the translational free energy contri- 
bution 18 . The vibrational free energy contribution will 
be of the order of the ZPE for the temperature range 
discussed here. As can be seen in Table U the ZPE value 
for CO2 is approximately 0.3 eV. Variations of AE mo1 in 
Eq. (|Sj) of this order of magnitude do not affect any of 
the conclusions discussed below, which justifies the rather 
crude approximation for the present purpose. 



C. Computational Setup 

With the approximations made in the last two sub- 
sections, the crucial quantities determining the sta- 
bility of surface structures are the total energies of 
the extended surfaces and of the involved gas phase 
molecules. These total energies are obtained by DFT 
calculations which have been performed within the full- 
potential (linearized) augmented plane wave + local or- 
bital (L)APW+lo metho d 21 ! 22 as implemented in the 
WIEN2k code 2 ^. 

All surfaces are simulated within the supercell ap- 
proach, using inversion symmetric slabs consisting of five 
Pd(100) layers with adsorption of oxygen and/or CO or 
the reconstructed surface oxide plus additional O/CO on 
both sides. The vacuum between consecutive slabs is 
at least 14 A. The adsorption layers and two outermost 
palladium layers have been fully relaxed. The muffin- 
tin radii have been set to it^T = 2.0bohr for palla- 
dium, Ry[T ~ l.Obohr for oxygen and Ry[ T = l.Obohr 
for carbon. Inside the muffin-tins the wave functions 
are expanded up to Z^ x = 12 and the potential up to 
Zg^ = 6. For the Pd(100)-(1 x 1) structure a [10 x 10 x 1] 
Monkhorst-Pack (MP) grid has been used to integrate 
the Brillouin zone (BZ). For larger surface unit cells the 
MP-grid has been reduced accordingly to assure an equiv- 
alent sampling of the BZ. Since calculations with dif- 
ferent MP-grids are not fully comparable, the energies 
of the clean metal surface, E^ (100) , and of bulk palla- 
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dium, needed to evaluate the Gibbs free energy of 

adsorption as defined in Eq. ([2|), have been calculated 
each time in the same supercell as the corresponding ad- 
layers or surface oxide structure. The energy cutoff for 
the expansion of the wave function in the interstitial is 
£™ f ax = 20 Ry and for the potential = 196 Ry. Us- 

ing these basis set parameters the average binding ener- 
gies per oxygen atom/ CO molecule are converged within 
50 meV, so that the Gibbs free energies of adsorption are 
converged within 1-5 me v/A 2 . 

All molecular calculations are done in rectangular su- 
percells with side lengths of (13 x 14 x 15)bohr (atoms) 
or (13 x 14 x 20) bohr (molecules) using only the T-point 
to sample the BZ. To obtain the binding energies of the 
three gas phase molecules, O2, CO and CO2, very accu- 
rately the cutoff for the expansion of the wave function in 
the interstitial has been increased to E^ ax = 37 Ry. With 
respect to this cutoff the binding energies are converged 
within 10-20 meV. 

The exchange-correlation (xc) energy is treated within 
the generalized gradient approximation (GGA) using the 
PBE 24 xc-functional. To assess the error introduced to 
the total energies by the approximate xc-energy, the most 
important calculations have also been repeated using an- 
other gradient corrected xc-functional, the RPBE 25 , as 
well as the local-density approximation, LDA2£. The in- 
fluence of these different xc-functionals on the results will 
be discussed below. 



III. O AND/OR CO ADSORPTION AT Pd(100) 

In this section we discuss the different adlayer struc- 
tures of oxygen and/or CO on Pd(100) that will be 
compared within the constrained thermodynamic equilib- 
rium approach. In addition to the on-surface adsorption 
of O/CO on Pd(100), the formation of a reconstructed 
(y/E x y/b)R27° (abbreviated with y/b in the following) 
surface oxide has experimentally been observed at higher 
oxygen exposur o 27 i 28 i 29 i 30 . Since the structure of the y/E 
surface oxide is quite different from the clean Pd(100) 
surface, also exhibiting different adsorption sites, we will 
discuss the two surfaces separately. 

A. Adlayers on Pd(100) 

1. Pure oxygen adsorption 

Experimentally, four different adsorption structures 
have been observed when exposing the Pd(100) surface 
to oxygen 2 ^ 2 ^^^!. These are a p(2 x 2), a c(2 x 2), 
a (5 x 5) and the y/E surface oxide structure. The y/E 
surface oxide will be discussed in detail in the next sec- 
tion. The (5 x 5) structure will not been considered in 
this study, since the structure and state of the oxygen 
atoms has not been well established so far, i.e. there ex- 
ists no well defined structural model. In addition, the 



(2y/2 X ^2)^45° (Sy/2 X V2)R45° (Ay/2 X y/2)BA&> 




FIG. 1: (Color online) Top views of the three experimentally 
characterized adlayer structures of CO on Pd(100). From left 
to right the coverage increases from O = 0.5 to O = 0.67 
and O = 0.75 ML. The yellow small spheres represent the CO 
molecules, the grey large spheres the Pd(100) surface (second 
layer atoms are darkened). 



(5 x 5) structure appears to be only of metastable char- 
acter and its formation was found to be very sensitive to 
the surface preparation and oxygen exposure range, so 
that in the oxidation as well as in the reduction process 
the (5 x 5) structure can be bypassed going directly from 
a (2 x 2) to a y/E structure and vice versa 30 . 

This leaves as ordered adlayers only the p(2 x 2) at 0.25 
monolayer (ML) coverage and the c(2 x 2) at 0.5 ML. 
Computing the binding energy of O atoms in all high- 
symmetry sites offered by the (100) surface (bridge, top 
and hollow), we find in both structures the fourfold hol- 
low site to be the most stable adsorption site. The bind- 
ing energies are calculated with respect to the O2 gas 
phase molecule (cf. Eq. (|HQ and the values of the bind- 
ing energies per oxygen atom are listed in Table HI1 It can 
be seen that the binding energy decreases going from the 
p(2 x 2) to the higher coverage c(2 x 2) adlayer, indicat- 
ing overall repulsive interactions between the adsorbed 
oxygen atoms. 



2. Pure CO adsorption 

The adsorption of CO on Pd(100) has been studied in- 
tensively^ 2 -^^^^^^! 3 -^. In all experimental studies 
it is found that CO binds upright in the bridge position 
via the C atom. At a coverage of O = 0.50 ML an ordered 
(2y/2 x y/2)R4b° adlayer is formed, which is compressed 
to a (3y/2 x y/2)R45° structure at 9 = 0.67 ML and a 
(Ay/2 x y/2)R45° at 6 = 0.75 ML (cf. Fig.HJ. The for- 
mation of the (2y/2 x y/2)R4b° adlayer rather than a 
simple c(2 x 2) structure is assumed to be mainly due 
to overall strongly repulsive interactions among the ad- 
sorbed CO molecules 40 . In the c(2 x 2) each CO molecule 
has four nearest neighbors at a distance of y/2a, where 
a is the length of the (lxl) surface unit cell. In the 
(2y/2 x v // 2)^45° structure, however, the CO molecules 
form a distorted hexagonal overlayer with only two near- 
est neighbors at a distance of y/2a and four neighbors 
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TABLE II: Calculated binding energies of on-surface adlayers 
of 0/CO on Pd(100). For the pure adlayer structures the 
binding energies are given per O atom resp. CO molecule. 
For the mixed structures the total binding energies are given. 
All values are in eV. 

Coverage 6 E hind 



(2x2) - top view 



Pure O structures: 

p(2 x 2)-O ho1 
c(2 x 2)-O ho1 

Pure CO structures: 

p(2 x 2)-CO br 

c(2 x 2)-CO br 

(2y/2 x v/2)i?45°-CO br 

(3^2 x V2)m5°-CO hr 

(4>/2 x V2)R4b°-CO hr 

(1 x l)-CO br 

Mixed O/CO structures: 

(2 x 2)-O hol -CO br 
(2 x 2)-20 hol -CO br 
(2 x 2)-O hol -2CO br 
(2 x 2)-20 hol -2CO br 



0.25 
0.50 



0.25 
0.50 
0.50 
0.67 
0.75 
1.00 



0.50 
0.75 
0.75 
1.00 



-1.35 
-1.10 



-1.92 
-1.92 
-1.92 
-1.73 
-1.63 
-1.31 



-2.95 
-2.93 
-3.96 
-3.64 



slightly further away at a distance of y 5 /2a. 

In addition to these experimentally observed structures 
also simple (lxl) and (2 x 2) adlayers with CO in top, 
bridge and hollow sites are considered in this work. The 
adsorption in bridge sites is found to be energetically fa- 
vored at all computed coverages. Table HT1 compiles there- 
fore only the results for bridge-bonded CO. Looking at 
the values in Table [III it can be seen that the binding 
energy per CO molecule is almost constant up to a cov- 
erage of O = 0.50 ML. For higher coverages the binding 
energy decreases due to overall repulsive interactions be- 
tween the adsorbed CO molecules. Comparing the two 
discussed structures with a coverage of 9 = 0.50 ML we 
in fact find them to be degenerate within our computa- 
tional accuracy (the (2V2 x V2)R45° is only favored by 
a few meV). These results are also fully consistent with 
previous DFT pseudopotential studies 25,40 . 



3. Mixed O/CO adsorption 

As to the adsorption of oxygen and CO on the Pd(100) 
surface, experimentally no ordered overlayers have been 
reported so far. If the Pd(100) surface is exposed to both 
oxygen and CO, the two adsorbed species tend to form 
separate domains instead 3 -. In their low energy electron 
diffraction (LEED) measurements Stuve et alM* observed 
that for a fully developed p(2 x 2)-O/Pd(100) structure 
the p(2 x 2) LEED pattern vanished after exposure to 
CO at a temperature of T = 80 K. Assuming a barrier of 
« 1.0 eV for the reaction of O and CO to form CO^^ it 
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FIG. 2: (Color online) Top views of the considered co- 
adsorbed overlayer structures of oxygen and CO on Pd(100). 
The small red (dark) spheres correspond to oxygen atoms, 
the small yellow (light) spheres represent the CO molecules 
and the large grey spheres the Pd(100) surface (second layer 
atoms are darkened). 



is rather unlikely that the adsorbed oxygen reacted with 
the CO at this low temperature. The disappearance of 
the LEED signal was thus interpreted as a CO induced 
disordering of the oxygen islands. To nevertheless obtain 
an idea about the simultaneous adsorption of O and CO 
on Pd(100) we set up different models for mixed overlayer 
structures that seemed to be most obvious from a com- 
binatorial point of view. For this we used (2 x 2) surface 
unit cells placing 1 to 2 O atoms and CO molecules in 
their favorite adsorption sites, i.e. O in hollow and CO in 
bridge sites. Excluding any structures where O and CO 
are closer to each other than the length a of the (lxl) 
surface unit cell we obtain the four different structures 
shown in Fig. [2j The corresponding average binding en- 
ergies calculated using Eq. ([3|) are listed in Table HI 

Similar to the pure adlayer structures of O or CO on 
Pd(100), the binding energy per adsorbate decreases with 
increasing coverage in the considered co-adsorption struc- 
tures. This trend in binding energies reflects the afore 
mentioned repulsive interactions among the adsorbed O 
atoms and CO molecules. For a coverage of O = 0.50 ML 
the binding energy of the mixed O/CO overlayer is less 
favorable by 70 meV than the sum of the binding energies 
of the respective c(2 x 2) pure adlayers of O and CO. This 
indicates that the repulsive interactions among adsorbed 
O and CO are even stronger than between O and O, re- 
spectively CO and CO in the pure adlayer structures. 
The adsorption of both O and CO in hollow sites in a 
c(2 x 2) structure further decreases the binding strength 
of the adsorbates. 

The even stronger repulsive interactions in these mixed 
adsorbate structures will thus favor the separation of oxy- 
gen and CO into separate domains. The adsorption of 
CO into ap(2x2)-0 structure can then induce a demixing 
of the two species and destroy the long-range order, which 
is consistent with the disappearance of the LEED spot 
intensities in the afore mentioned experiment by Stuve et 
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surface oxide 



FIG. 3: (Color online) Top view of the (y/E x y/E)R27° surface 
oxide structure (left part) and the Pd(100) substrate (right 
part) and their corresponding surface unit cells (black lines). 
Red small spheres represent the lattice oxygen atoms, large 
light-blue ones the reconstructed Pd atoms and large dark- 
blue spheres Pd atoms belong to the underlying Pd(100) sub- 
strate. The arrow indicates the easy shift direction of the 
surface oxide trilayer along the Pd(100) substrate. 



B. O and CO on the surface oxide 

At a coverage of 6 = 0.80 ML, the adsorbed oxygen 
induces a reconstruction of the Pd(100) surface forming 
the y/E surface oxide structure. An atomistic model for 
the surface oxide (cf. Fig. [3]) was proposed in a detailed 
experimental and theoretical study, revealing that the 
structure can actually be described as a PdO(lOl) tri- 
layer on Pd(100) 4 -. The surface unit cell contains four 
palladium and four oxygen atoms. Two Pd atoms are 
fourfold and the other two are twofold coordinated by 
oxygen. There are also two kinds of O atoms. Two O 
atoms sit on-top of the reconstructed Pd layer and two at 
the interface to the Pd(100) substrate forming the afore 
mentioned trilayer structure. 

In addition to the previous discussion on this surface 
oxide structure^, we note that the potential energy sur- 
face for a registry shift of the entire surface oxide trilayer 
over the Pd(100) substrate is very shallow. In particular, 
shifts of up to ~ 0.5 A along the direction as indicated by 
the arrow in Fig. [3] and parallel to the surface lead only to 
energy variations of less than 50 meV per y/E surface unit 
cell. Without anchoring by e.g. defects, the lateral posi- 
tion of the ideal surface oxide over the Pd(100) surface is 
thus not well denned. In extended test calculations, we 
verified that this uncertainty in the lateral position has 
no consequences on the adsorption energetics discussed 
here and therefore simply use in the present work the 
lateral position determined in Ref. [43|. 

For the additional adsorption of oxygen and/or CO on 
the surface oxide structure there are several possibilities. 
As can be seen in Fig. [H the surface oxide exhibits top, 
bridge and hollow high-symmetry sites. Due to the sym- 
metry of the underlying Pd(100) substrate multiple ad- 
sorption sites of the same type (top, bridge or hollow) are 
not fully equivalent, but still very similar. Only the top 
sites differ substantially depending on whether the cor- 
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FIG. 4: (Color online) Top views of high- symmetry adsorp- 
tion sites on the (y/E x y/E)R27° surface oxide structure. Red 
small spheres represent the lattice oxygen atoms, large light- 
blue ones the reconstructed Pd atoms and large dark-blue 
spheres Pd atoms belong to the underlying Pd(100) substrate. 



TABLE III: Calculated binding energies of CO on the y/E 
surface oxide structure calculated using Eq. ©. 2CO denotes 
that two CO molecules are adsorbed in the y/E surface unit 
cell. The adsorption of two CO molecules in top4f sites or one 
CO in bridge and a second in a top2f site does not lead to a 
stable structure. All values are per CO molecule and in eV. 





E-ibind 




E-ibind 




E-ibind 


C Qbr 

QQtop2f 

C Qtop4f 


-0.93 
-0.62 
-0.13 


2CO br 

2C Qto P 2f 
2C Qto P 4f 


-0.75 
-0.51 


2QQbr,top2f 
2C Qbr,top4f 
2£JQtop2f ,top4f 


-0.45 
-0.29 



responding Pd atom is twofold or fourfold coordinated 
by oxygen. As yet there is no experimental information 
available for the additional adsorption of O or CO on 
the y/E surface oxide. Therefore, we performed a sys- 
tematic investigation of possible overlayer structures of 
O and CO involving the reconstructed y/E surface oxide 
and starting with the adsorption sites depicted in Fig. [H 
In a first step only one CO molecule is adsorbed in 
any of the ten adsorption sites. We find that the four 
different hollow sites are not stable, i.e. upon relax- 
ation the adsorbed CO molecule always moves to the 
respective neighboring bridge site. The binding energies 
of the two bridge (br I and br II) sites are almost de- 
generate as expected from the only slight structural dif- 
ferences, i.e. the binding energies as defined in Eq. ([5]) 
differ by less than O.leV. Similarly, the two top sites at 
the two twofold (top2f I and top2f II) oxygen coordi- 
nated Pd atoms exhibit nearly equivalent binding ener- 
gies (Ai^ bmd < 50meV). The same is also found for the 
two top sites at the two fourfold (top4f I and top4f II) 
oxygen coordinated Pd atoms. In the following we will 
thus treat these correspondingly very similar sites as be- 
ing identical. The bridge site is the most stable adsorp- 
tion site followed by the top2f site, and on the top4f site 
the CO molecule is already only weakly bound (cf. Ta- 
ble [Till). 

It is also possible to adsorb two CO molecules in the 
y/E surface unit cell, and again we find the adsorption in 
the two bridge sites most stable. Mixing the adsorption 
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sites, i.e. placing e.g. one CO in a bridge site and a 
second one in a top2f site does always yield less stable 
structures than the adsorption in like sites (cf. Table HTT|). 
It also becomes energetically very unfavorable to place 
more than two CO molecules in the y/b surface unit cell as 
deduced from test calculations with up to 4 CO molecules 
per surface unit cell. 

The adsorption of oxygen in the high-symmetry sites 
of the y/b surface oxide gives similar results. Again, the 
bridge site is the most stable site (^qSvs = — 0-14 eV) 
and the hollow sites are unstable upon relaxation. How- 
ever, it is now not possible to adsorb any oxygen in the 
top2f and top4f sites in the y/b unit cell. If one CO 
molecule and one O atom are adsorbed simultaneously, 
likewise the adsorption in two bridge sites is found to be 
favored. 

Our calculations with one and two adsorbates per y/b 
surface unit cell thus indicate noticeable overall repulsive 
lateral interactions. We correspondingly extended our 
calculations also to more sparse adlayers in (2 x 1) and 
(1x2) y/b surface unit cells, but did not find significant 
lateral interactions extending across the y/E surface unit 
cell. The obtained binding energies were always to within 
25meV per adsorbate of those calculated for the same 
adsorption site in a (1 x 1) cell. 

To increase the configuration space of considered over- 
layers based on the y/E surface oxide, we also considered 
structures, in which the original surface oxide structure 
is slightly modified. In a first step, one of the upper 
hollow site oxygen atoms is removed, and CO and O 
are again placed in the afore described sites. Also in 
this modified structure the bridge site results as the pre- 
ferred adsorption site. If both upper oxygen atoms are 
removed from the y/E surface oxide, we find the struc- 
ture already to some extend destabilized. Due to the 
lower palladium density in the reconstructed layer (4 Pd 
atoms on 5 Pd(100) substrate atoms per unit cell) the 
structure is rather open, and the structural relaxation 
showed that the Pd atoms can now move quite easily in 
lateral direction on the surface. It is still possible to ad- 
sorb O/CO in any of the other adsorption sites, but we 
found these structures either unstable upon relaxation or 
in general much less stable than the corresponding struc- 
tures, where the original O hollow vacancy is refilled. 

In addition to simply removing the upper oxygen 
atoms from the y/E surface oxide structure and filling ad- 
ditional on-surface adsorption sites, we also substituted 
the topmost O atoms by CO molecules, and placed ad- 
ditional O atoms and/or CO molecules in top2f, top4f 
and bridge sites. In the substituted hollow site the CO 
binds quite strongly. However, since the O atoms bind 
even stronger to this site by 0.6 eV, O will always prefer- 
entially occupy the hollow sites when the two adsorbates 
compete for this sites. The situation is reversed for ad- 
sorption in the bridge sites. Regardless of whether O or 
CO occupies the threefold hollow sites, we always find a 
stronger binding by at least 0.7 eV of CO compared to O 
at these bridge sites. 



Even if we leave out the four unstable hollow sites in 
Fig. 2] there is finally still a huge amount of possible over- 
layer structures that can be created by combinatorially 
placing an arbitrary number of O and CO per surface 
unit cell into any of the available sites. Since there is 
only little known about the adsorption on the surface 
oxide none of these structures can a priori be excluded. 
Motivated by the strong repulsive interactions seen in 
our calculations with one or two adsorbates per surface 
unit cell, we nevertheless discard quite a number of these 
structures with the criterion that no two adsorbates may 
sit in directly neighboring sites (i.e. at a distance of less 
than 1/4 of the length of y/E surface unit cell). This 
still leaves 92 "plausible" structures and DFT calcula- 
tions were performed for all of them (for more details on 
the calculated overlayer structures see Ref. 44). Out of 
these only 55 are stable upon relaxation and are then 
considered in our constrained "phase diagram" . 



IV. STABILITY IN A CONSTRAINED 
EQUILIBRIUM WITH AN 2 AND CO GAS 
PHASE 

Using the Gibbs free energy of adsorption as defined in 
Eq. ([2]) it is possible to compare the stability of all calcu- 
lated adsorption structures in a constrained equilibrium 
with an O2 and CO gas phase. In this approach the sur- 
face is considered to be in equilibrium with two separate 
gas phase reservoirs of O2 and CO, i.e. the formation 
of CO2 is neglected in the gas phase as well as on the 
surface 16 . In the gas phase this assumption seems quite 
justified, since the direct reaction of O2 and CO is ki- 
netically hindered by a huge free energy barrier. On the 
surface, though, the reaction is actually supposed to take 
place at the working catalyst. Here, the assumption of 
a constrained equilibrium is thus only reasonable as long 
as the O2 and CO adsorption and desorption events are 
much more frequent than the reaction, so that the sur- 
face can maintain its equilibrium with the two gas phase 
components. In other words, the approach amounts to 
assuming that the kinetics of the on-going catalytic reac- 
tions does not change the surface structure and composi- 
tion. The resulting constrained surface "phase diagram" 
can thus only provide a first idea of the possibly relevant 
surface structures, which needs to be scrutinized explic- 
itly treating the surface kinetics. On the other hand the 
constrained equilibrium approach allows to rapidly screen 
a huge number of (structurally quite different) surface 
configurations and compare their stability over a wide 
(T,p)-range of possible gas phase conditions. By using 
ab initio thermodynamics we are thus able to identify 
those regions in (T,p)-space that are likely catalytically 
active and where we then zoom in with more sophisti- 
cated methods. This second, refining step explicitly in- 
cludes the effect of the reaction kinetics on the average 
surface composition (using kinetic Monte Carlo simula- 
tions) and the results of these simulations, which are re- 
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stricted to the most relevant surface structural models 
identified here, are discussed in a second paper. This 
will then also bring us into a position to scrutinize the 
validity of the assumption of the constrained equilibrium 
approach employed here. 



A. Surface Phase Diagram 

The results of the above describe calculations concern- 
ing the afore discussed over layers of oxygen and CO on 
Pd(100), also including the (modified) y/E surface ox- 
ide structure, are summarized in Fig. which shows the 
most stable "phases" (i.e. the ones that maximize AG ads 
as defined in Eq. ©) for any given chemical potential of 
the oxygen and CO gas phase. The dependence on the 
chemical potentials has also been converted into more 
intuitive pressures scales for temperatures of T = 300 K 
and T = 600 K in the top two x-axes for oxygen and 
the right two ?/-axes for CO. As expected from the 
Gibbs phase rule, stable "phases" cover volumes in the 
3-dimensional space spanned by (T,po 2 ?Pco) 7 and coex- 
istence regions between two "phases" cover areas. Cor- 
respondingly in the 2-dimensional sectional plane shown 
in Fig. [5] this leads to areas and lines for the two cases, 
respectively. Out of the set of 191 tested structures only 
11 turn out to be "phases" appearing in Fig. [5j 

We will start the discussion of the obtained constrained 
surface "phase diagram" in Fig.[5]in the lower left corner. 
Here, both the chemical potential of oxygen and CO are 
very low, i.e. the oxygen/CO content in the gas phase is 
insignificant and consequently the clean Pd(100) surface 
results as the most stable system state. If we now move 
along the x-axis in Fig. [5] to the right we reach more and 
more oxygen-rich conditions, while the CO content in the 
gas phase is kept low. The increase in the oxygen chemi- 
cal potential leads to a stabilization of oxygen containing 
structures with increasing coverage: First, the p(2 x 2)-0 
over layer on Pd(100) 45 , then the y/E surface oxide and 
finally the PdO bulk oxide representing thick bulk-like 
oxide films on the surface. This sequence of stable struc- 
tures was also confirmed by in situ surface X-ray diffrac- 
tion measurements 4 . Interestingly, the c(2 x 2) structure 
observed under UHV condition s 27 ! 28 ! 29 ! 30 * 3 * 1 does not ap- 
pear in the surface "phase diagram" . Thus we conclude 
that the c(2 x 2) structure is most likely a meta-stable 
state produced by the exposure kinetics. 

If the oxygen content is kept low in the gas phase and 
the CO content is gradually increased, i.e. moving from 
the lower left corner along the y-axis to the top in Fig. 
a series of ordered CO adlayer structures with increasing 
coverage is stabilized on Pd(100). At first we find the 
also experimentally observed (2y/2 x y/~2)R4:5 , (3y/2 x 
y/2)R4b° and (4\/2 x y/2)R4b° structures, and finally for 
a very high CO content in the gas phase a (1 x 1) structure 
with 1 ML CO in bridge sites. 

Starting again in the lower left corner of Fig. [5] we now 
move along the diagonal, which corresponds to increas- 



ing both the oxygen and CO content in the gas phase. 
Intuitively, one would expect co-adsorption structures of 
O and CO on Pd(100) to become favorable. But none 
of the above discussed ordered co-adsorption structures 
(cf. Fig. [2]) are found to be a most stable "phase" under 
any gas phase conditions. This is consistent with the al- 
ready mentioned experimental findings, that O and CO 
prefer to form separate domains rather than ordered co- 
adsorbed overlayers 31 . The lower stability of such mixed 
structures can be explained by the strongly repulsive in- 
teractions between adsorbed O and CO, which lead to 
a significant decrease in binding energies (cf. Table HI]) . 
However, we can of course not exclude that ordered ar- 
rangements with different periodicities than those con- 
sidered here would not lead to a lowering in the repulsive 
interactions. To take a reasonable part in the "phase di- 
agram" , though, the binding energies would have to in- 
crease by as much as 0.3-0.5 eV per O atom/ CO molecule 
compared to the now proposed structures. Compared 
to the pure overlayer structures this would even imply 
the necessity of attractive interactions between the two 
species. 

For high oxygen and CO content in the gas phase (up- 
per right part of Fig. [5]) we find instead co-adsorbed 
structures involving the y/E surface oxide (two orange- 
white hatched regions) to become stable. These two 
mixed structures correspond to the surface oxide with 
one and two CO molecules adsorbed in bridge sites, re- 
spectively. In a small range of very oxygen-rich and in- 
termediate CO gas phase conditions additional O adsorp- 
tion in the bridge sites of the y/E structure leads also to 
a stable "phase" (dark red- white hatched region). 

Looking again at the whole surface "phase diagram" in 
Fig. [5] the 11 "phases" can be divided into three groups: 
on-surface adlayer structures of O or CO on Pd(100), 
phases involving the y/E surface oxide structure, and the 
stability region of the bulk oxide as discussed in sec- 
tion III Bt Focusing on these three different groups of 
phases two important conclusions with respect to the rel- 
evance of oxide formation under catalytic reaction con- 
ditions can be drawn. First, the formation of a thick, 
bulk-like oxide (grey cross-hatched area) at the surface 
under technologically relevant gas phase conditions of O2 
and CO (black bar in Fig. El pi ~ 1 atm, T ~ 300- 600 K) 
can be ruled out. This is thus in marked contrast to the 
much wider stability range of bulk Ru02, which does ex- 
tend to these conditions 18 . Second, the stability region of 
the y/E surface oxide (hatched area) does extend to such 
conditions. In fact, it is either this monolayer thin surface 
oxide or CO adlayers on Pd(100), which are neighboring 
"phases" around the catalytically active region in (T,_p)- 
space. As can be seen in Fig. [5] the catalytically relevant 
conditions are actually right at the boundary between the 
y/E surface oxide structures and the CO overlayer struc- 
tures on Pd(100). Small changes in the Gibbs free en- 
ergy of adsorption causing a shift in this boundary could 
thus well affect the conclusion as to which is the most 
stable "phase" under these conditions. Apart from the 
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FIG. 5: (Color online) Surface "phase diagram" of the Pd(100) surface in constrained thermodynamic equilibrium with an O2 
and CO gas phase. The atomic structures underlying the various stable (co-) adsorption structures on the metal and the surface 
oxide, as well as a thick bulk-like oxide film (indicated by the bulk unit cell) are also shown (Pd = large blue spheres, O = 
small red spheres and C = small yellow spheres). The black bar marks gas phase conditions representative of technological CO 
oxidation catalysis, i.e. partial pressures of 1 atm and temperatures of 300 — 600 K. The dependence of the chemical potentials 
of the two gas phases is translated into pressure scales for T = 300 K and T = 600 K (upper x-axes and right 2/- axes). 



assumption of a constrained thermodynamic equilibrium, 
the main uncertainties, which may cause such changes in 
AG ads are the approximate DFT total energies and the 
neglected free energy contributions. We verified that the 
uncertainties in the DFT total energies due to the numer- 
ical approximations (supercell setup, finite basis set) are 
not significant in this respect. What needs to be scruti- 
nized are therefore the uncertainties due to the approx- 
imate xc-functional underlying the DFT total energies 
and the neglected free energy contributions to AG ads . 
This will be done in the next two subsections. 



B. Evaluating the Gibbs Free Energy 

To calculate the Gibbs free energy of adsorption AG ads 
the Gibbs free energies of the different components enter- 
ing Eq. ([2]) have to be evaluated. The Gibbs free energy 

G(T, p) = E tot + F vib - TS conf + pV (9) 

comprises contributions from the total energy E tot , the 
vibrational free energy F vlb including the zero point en- 
ergy (ZPE), the configurational entropy g conf and the 
pV-term. In Eq. ([2]) we substituted the Gibbs free en- 
ergy difference by only the leading term, the total energy 
differences, which can be directly obtained from DFT 
calculations. We assess the uncertainty introduced by 
this approximation by an order of magnitude estimate 
of the remaining contributions to AG ads . If this first 
approximation reveals that the results are significantly 



influenced by considering all contributions to the Gibbs 
free energy, the respective terms have to be calculated ex- 
plicitly. However, the order of magnitude estimate can be 
obtained very easily and is thus helpful to decide whether 
or not it becomes necessary to evaluate the entire Gibbs 
free energy. 

Following the discussion in Refs. [l3| and [l6| the con- 
tributions to the Gibbs free energy of adsorption arising 
from the pV-term and the configurational entropy are 
rather small for systems like the one presented here. We 
will thus only discuss the most crucial approximation, 
namely the neglect of the vibrational free energy contri- 
bution, AF ads ' vlb to the Gibbs free energy of adsorption. 
In order to estimate the size of AF ads ' vlb we approximate 
the phonon density of states (DOS) within the Einstein 
model by one characteristic frequency. Following the ap- 
proach outlined in Ref. [l3| the vibrational contribution 
to the Gibbs free energy of adsorption for the p(2 x 2)- 
O/Pd(100) structure can then be estimated as 

AF ads ' vib (T) * (10) 

~ ^ I Z^vib/^i -,surf \ 1 TTwib,ZPE/ -,gas\ \ 

where F vib (T, u) is the frequency dependent function 

F vib (7» = l&j + k B T\n(l - e~P huj ) (11) 
with (3 = l/k B T. 

^o 1 — Pd * s characteristic vibrational 
frequency of an oxygen atom adsorbed in a fourfold hol- 
low site on Pd(100), and Co^ is the stretch frequency 
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contribution are included. There are some small shifts 
in the boundaries between stable phases, but none of the 
stable structures disappears from and none of the unsta- 
ble ones appears in the "phase diagram". Furthermore, 
the boundary between the surface oxide structures and 
the CO adlayers on Pd(100) is only marginally affected 
and consequently remains in the vicinity of the catalytic 
active region in (T,p) -space. We are thus confident that 
for the here discussed results the approximation of the 
Gibbs free energy differences by only the total energy 
terms in Eq. ([2]) is justified. 



FIG. 6: (Color online) Vibrational contribution to the Gibbs 
free energy of adsorption for the p(2 x 2)-0 adlayer on 
Pd(100), within the Einstein approximation using different 
characteristic frequencies (see "inset" and text). 



of the O2 gas phase molecule. The change in the vibra- 
tional contribution of the Pd atoms in the clean and oxy- 
gen covered Pd(100) is neglected, so that the vibrational 
contribution to AG ads is given by the difference in the 
vibrational energy of the O2 molecule in the gas phase 
and the oxygen atom on the surface. For the oxygen 
molecule only the ZPE (not included in the presented 
DFT total energies) has to be considered, since all re- 
maining contributions to F vlb are already contained in 
A/io (cf. Eq. ([2])). In Fig. [6] the resulting vibrational 
contribution AF ads ' vlb (T) is shown for a characteristic 
frequency of ^q^im = 48meV for the Pd-0 stretch fre- 
quency of an adsorbed O atorn^ and an 0-0 vibrational 
frequency of Cu^ = 196 me V— . We find that for the 
p(2 x 2)-O/Pd(100) structure the vibrational contribu- 
tion to AG ads stays below ~ 3meV/A 2 for temperatures 
up to T = 1000 K (black-solid line in Fig. [6]). Since the 
chosen c^o-Pd ^ s oru y a rou gh guess, we also allow for a 
±50 % change of this value, but even then the contribu- 
tion does not increase considerably (red-dashed, —50% 
and blue-dotted, +50%, lines in Fig. [6j). 

Similar estimates have been performed for every struc- 
ture considered for the surface "phase diagram" in Fig.[5j 
For CO containing structures AF ads ' vlb comprises two 
different contributions when comparing the gas phase 
and the adsorbed state, namely the change of the C- 
O stretch vibration due to the adsorption and an addi- 
tional Pd-C vibration. For structures involving the y/b 
surface oxide an additional contribution arises from the 
change in the vibrational energy between bulk and sur- 
face Pd atoms. The y/b structure contains one surface 
atom less than the corresponding Pd(100) surface per 
unit cell which also has to be balanced by the bulk reser- 
voir (A7Vp d = -1 in Eq. ©). 

With this procedure we find that for all structures 
AF ads ' vlb stays always below 10meV/A 2 for tempera- 
tures up to T = 600 K. The surface "phase diagram" dis- 
cussed in the previous section is not significantly changed, 
if these estimated, maximum values for the vibrational 



C. Influence of the xc-functional 

To obtain an estimate of the uncertainty in Fig. [5] due 
to the choice of the PBE-GGA xc-functional we reeval- 
uate the constrained surface "phase diagram" using the 
RPBE 25 and LDA2£ xc-functionals. The corresponding 
surface "phase diagrams" are shown in Fig. and are to 
be compared with Fig. [5l The "phase diagram" obtained 
by using the RPBE approximation for the xc-energy (top 
graph, Fig. Ufa)) looks in fact very similar to the previ- 
ously discussed one. There are some shifts in the ac- 
tual phase boundaries, and the stability regions of the 
high coverage phases of oxygen and CO, the y/b-\-O hr 
and (1 x l)-CO br /Pd(100) structures, are shifted outside 
the shown range of chemical potentials. Yet, the overall 
topology is fully conserved. In agreement with the PBE 
results the bulk oxide is not a stable "phase" anywhere 
near ambient gas phase conditions. Due to the smaller 
heat of formation (cf. Table. [Ij), the stability region of the 
bulk oxide is in fact much further away from this region. 
Most importantly, the boundary between the simple over- 
layer structures on Pd(100) and the phases based on the 
y/b surface oxide structure (hatched area) is only very 
little influenced, so that the most interesting region for 
oxidation catalysis (black bar) is again right at the tran- 
sition between a CO covered Pd(100) surface and phases 
involving the surface oxide. 

Even if the LDA is used as xc-functional (bottom plot, 
Fig. 0(b)) the stability region of the bulk oxide does not 
extend to the catalytically relevant gas phase conditions. 
It covers a much larger range compared to the "phase di- 
agrams" obtained with the two GGA functionals though. 
It should also be noted that in the shown LDA "phase di- 
agram" the range of O and CO chemical potentials is en- 
larged to lower values to include the stability region of the 
clean metal surface, which in the LDA appears at much 
lower gas phase concentrations. This is a consequence 
of the strong overbinding in the LDA, which stabilizes 
any adsorbate structure at much lower pressures at the 
surface. Nevertheless, the range of technologically rele- 
vant catalytic gas phase conditions (black bar) lies again 
right at the boundary between the CO covered Pd(100) 
surface and y/b surface oxide structures. 

From these results, we conclude that the choice of the 
xc-functional does have a strong influence on the abso- 
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FIG. 7: (Color online) Surface "phase diagram" for the 
Pd(100) surface in a constrained equilibrium with an O2 and 
CO gas phase. The top plot shows the RPBE results and the 
bottom one the LDA results. The color coding of the different 
phases is equivalent to the one used in Fig. [5] The black bar 
marks again gas phase conditions representative of technolog- 
ical CO oxidation catalysis, i.e. partial pressures of 1 atm and 
temperatures between 300-600 K. Note that the shown range 
of chemical potentials is larger for the LDA results, compared 
to the RPBE and the range in Fig. [5] 



lute values of the binding energies, and thereby also on 
the location of most of the phase boundaries in Fig. [5] 
and [3 However, particularly the boundary between the 
CO over layer structures on Pd(100) and the y/b surface 
oxide is almost unchanged. The conclusion on the prox- 
imity of both phases to the catalytically relevant gas 
phase conditions seems therefore untouched by the un- 
certainty due to the approximate xc-functional. 



V. CONCLUSIONS 

The stability of the Pd(100) surface has been investi- 
gated in a constrained thermodynamic equilibrium with 
a two component gas phase consisting of O2 and CO. In 
this approach the formation of CO2 in the gas phase and 
at the surface is not considered, such that the effect of the 
surrounding gas phase on the surface structure and com- 
position is modeled (as a first approximation) through 
the contact with independent reservoirs representing the 



reactants. To establish the actual surface "phase dia- 
gram" a large number of different structures with O and 
CO adsorbed in high symmetry sites on the Pd(100) sur- 
face and on the y/b surface oxide structure have been 
considered. We find that under gas phase conditions of 
ambient temperatures and pressures, as applied in het- 
erogenous oxidation catalysis, it is either the nanometer 
thin surface oxide structure or a CO covered Pd(100) sur- 
face that is stable, whereas the stability region of the bulk 
oxide does not extend to these gas phase conditions 18 . 

To obtain an estimate of the uncertainty introduced by 
the choice of the specific xc-functional the total energies 
entering the surface "phase diagram" have been calcu- 
lated using the PBE, RPBE and LDA. Comparing the 
"phase diagrams" for these three different xc-functionals, 
partly dramatic shifts in the positions of the boundaries 
between the different stable phases can be observed. Yet, 
it can also be seen that the position of the boundary be- 
tween the CO covered Pd(100) phases and the y/b surface 
oxide structures is in fact little affected. For all three xc- 
functionals, the catalytically active region is very close 
to this boundary, but just still within the stability range 
of the surface oxide. We also verified that this finding is 
not affected by the numerical uncertainties and the free 
energy contributions neglected in our approach. 

These results suggest that the only monolayer thin sur- 
face oxide structure might indeed be a relevant structure 
for the CO oxidation reaction on Pd(100) at technological 
relevant pressures. It has to be considered, though, that 
there are still two notable approximations in the con- 
strained atomistic thermodynamics approach, as it has 
been applied here. First, configurational entropy is not 
included, and at finite temperatures a first effect would 
be to smear out the phase boundaries in the "phase di- 
agram" and create coexistence regions, where the ener- 
getically low lying configurations are mixed according to 
the law of mass action. At the phase boundary relevant 
for CO oxidation catalysis, this could e.g. lead to the 
formation of coexisting domains of CO covered Pd(100) 
and surface oxide patches. Second, the neglected kinet- 
ics of the on-going catalytic CO2 formation might signif- 
icantly change the stability range of the different phases 
as obtained in the constrained equilibrium approach. A 
possible coexistence between patches of metallic Pd(100) 
covered by the reactants, and patches of the surface oxide 
(possibly also with adsorbed CO) could then even under 
steady-state conditions go hand in hand with a contin- 
uous formation and decomposition of the oxidic phase. 
Such oscillations in the morphology of the catalyst sur- 
face could again significantly influence the catalytic func- 
tion of this surface. 

Based on the insight described in this paper we thus 
identify either the y/b surface oxide or CO covered 
Pd(100) as surface structures most relevant under cat- 
alytically interesting gas phase conditions. It is therefore 
specifically the stability of the y/b surface oxide against 
CO induced decomposition, which needs to be scrutinized 
to assess the active state of the surface under reaction 
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conditions. In a second paper following the present one, 
we correspondingly extend our study to investigate this 
stability. Using kinetic Monte Carlo simulations we then 
explicitly scrutinize the two major uncertainties of the 
present approach, namely the effect of configurational 
disorder and the kinetic effects due to the on-going reac- 
tions. 
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